Astronomy & Astrophysics manuscript no. final 


February 4, 2008 


(DOI: will be inserted by hand later) 





o 
o 

(N 
>^ 



> 

oo 

vn 
o 
\o 
o 

o 



Gamma-ray binaries: pulsars in disguise ? 

Guillaume Dubus^'^ 

Laboratoire Leprince-Ringuet, UMR 7638 CNRS, Ecole Poly technique, 91128 Palaiseau, France; e-mail: dubus(ain2p3 . fr 
Institut d' Astrophysique de Paris, UMR 7095 CNRS, Universite Pierre & Marie Curie, Paris 6, 98 bis bd Arago, 75014 Paris, 
France 



Draft Februai-y 4, 2008 

Abstract. LS 5039 and LSI +61° 303 are unique amongst high-mass X-ray binaries (HMXB) for their spatially-resolved radio 
emission and their counterpart at >GeV gamma-ray energies, canonically attributed to non-thermal particles in an accretion- 
powered relativistic jet. The only other HMXB known to emit very high energy (VHE) gamma-rays, PSR B 1259-63, harbours 
a non-accreting millisecond pulsar. 

The purpose is to investigate whether the interaction of the relativistic wind from a young pulsar with the wind from its stel- 
lar companion, as in PSR B 1259-63, constitutes a viable scenario to explain the observations of LS 5039 and LSI +61°303. 
Emission arises from the shocked pulsar wind material, which then flows away to large distances in a comet-shape tail, repro- 
ducing on a smaller scale what is observed in isolated, high motion pulsars interacting with the interstellar medium. 
The timescales for acceleration and radiation of particles at the shock between the pulsar wind and stellar wind are calculated. 
Simple expectations for the spectral energy distribution (SED) are derived and are shown to depend on few input parameters. 
Detailed modelling of the particle evolution is attempted and compared to the observations from radio to TeV energies. 
Acceleration at the shock provides high energy electrons that steadily emit synchrotron in X-rays and inverse Compton scatter 
stellar light to y-rays. Electrons streaming out of the system emit at IR frequencies and below. The overall aspect of the SEDs 
is adequately reproduced for standard values of the parameters. The morphology of the radio tail can mimic a microquasar 
jet. Good agreement is found with the published VLBI map of LS 5039 and predictions are made on the expected change in 
appearance with orbital phase. 

The pulsar wind scenario provides a common, viable framework to interpret the emission from all three y-ray binaries. 

Key words, acceleration of particles — binaries: close — pulsars: general — ISM: jets and outflows — Gamma rays: theory 
— X-rays: binaries 

1 . Introduction luminal ejections can be taken as different manifestations of a 

unique, underlying ejection mechanism. 

X-ray binaries are systems composed of a neutron star or black 

hole in orbit around a high or low-mass normal star (HMXBs Similar types of radio emission (compact or superluminal 

and LMXBs). Non-thermal radio emission has been detected ejections) are seen in active galactic nuclei (AGN), prompt- 

in X-ra y binaries and is usually associated with outflows (see ing speculations that phenomenology known to one type of 

iFendeii p006, for a review). In a few objects, radio observa- object may have its counterpart in the other. The detection of 

tions are able to follow the propagation of radio emission up ^ery high energy y-rays from the X-ray binary LS 5039 by the 

to parsec scales, after initial flaring in the central source. The HESS collaboral.ion would appear to vindicate this conjecture 

inferred speeds of the collimated outflow are relativistic with ( Aharonian et aU|2005a|). Observations in the 1990s have es- 

bulk Lorentz factors F of a few. These apparently superluminal tabhshed that blazars, a subtype of AGN, emit high energy (HE 

ejection events have been tentatively a ssociated with trans itions GeV) and very high energy (VHE, TeV) y-rays (jHoffman etalj 
from a hard to a soft X-ray spectrum JFender et allEooi . Blazars have their jet aligned close to the line-of-sight, 

, , . n resulting in a very strong relativistic boosting of the iet non- 

In many more objects, the radio emission has a flat spec- 
trum and is correlated with a hard X-ray spectral state. Usually 
the radio emission is unresolved but imaging by Very Long 

BaseUne Interferometry (VLBI) of Cyg X-1 in this state has re- LS 5039 is composed of a 06.5V star and an unidenti 



thermal emission. The same configuration could be found in an 
X-ray binary and this is perhaps the case in LS 5039. 



veal ed a compact, colUni ated radio outflow on a scale of a few fled compa c t object in a 3.9 day or bit (Motc h etall 1 1997 



AU jStirling etalll200lh . Asymmetries in the surface bright- Clark et al., 2001; McSwain & GiesL [20021: fMcSwain et al 



ness probably reflect moderate relativistic boosting (v/c ~ 0.3). .2004:.Casares et al.. ■2005 b). Much of the interest in this object 
It is as yet unclear whether compact radio emission and super- stems from its tentative association with the EGRET source 
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3EG J 1824- 15 14 and from the fact that its radi o emission was 
shown to be extended on a milliarcsecond scale dParedes et all 
l200d ). By analogy with blazars, it is tempting to attribute the 
radio and y-ray emission to high energy particles accelerated in 
a jet. LS 5039 is very similar to LSI H-61°303 an X-ray binary 
with a BOVe star in a 26.5 day orbit whose association with a 
high energy source dates back to COS-B (Gregory & Taylor, 
iHntchings & Cramntonl \mA ICasares et all l2005all. 
The associations rely purely on positional coincidence and 
this is fairly poor, the EGRET 99% confidence levels being 
roughly at 0.5° for LS 5039 and 0.2° for LSI -i-61°303 (3EG 
J0241H-6103). The HESS detection of VHE emission coinci- 
dent to within an arcmin with the position of LS 5039 identi- 
fies this X-ray binary as the source of y-rays. By extension, it 
places the EGRET tentative identification on a firmer footing. 

VHE y-ray emission is at least evidence that particles 
can be accelerated to very high energies in X-ray binaries. 
However, is it necessarily associated with blazar-like model, 
involving accretion of mate rial, ejection and non- therm al accel- 
eration in a relativistic iet? lMaraschi & TrevesI ( Il98ll) pointed 
out early on that the spindown of a young, rapidly-rotating pul- 
sar could power the y-ray emission of LSI +61°303. The rela- 
tivistic pulsar wind would be confined by the stellar wind from 
the massive companion. Particles accelerated at the termination 
shock then produce the non-thermal emission. 

There is at least one system where such a scenario is un- 
doubtedly at work: PSR B1259-63 (=SS 2883), a 47.7 mil- 
lisec ond (ms) radio pulsar in a 3.4 yr orbit aroun d a B2Ve 
star ([.Johnston et all 1 19941: [Manchester et all Ll995). Variable 
y-ray emission was detected by HESS during periastron pas- 
sage (l^aronian et al., 2005b). In LS 5039 and LSI 4-61°303, 
the discovery of resolved radio emission was interpreted as 
the signature of a relativistic jet and recently, the attention 
seems to have focused on accretion/e jection scenarios (but see 
iLea hv 2004; ' Martocchia et Zl l2004 . The pulsar scenario has 
become unwonted. Yet, pulsars interacting with their surround- 
ing medium can also lead to large scale emission, the bow 
shock nebula of high velocity m s pulsars even taking a comet 
tail appearance (e.g. iGaensleJ. f2005 ). The resolved 'jets' of 
LS 5039 and LSI -i-61°303 might be small-scale analogues, but 
with the trail originating from the interaction with a companion 
wind rather than with the ISM. 

The similarities between LS 5039 and LSI +61°303 have 
long been noticed: resolved radio emission, hard steady X- 
ray fluxes and EGRET emission in the GeV range at a level 
~ 10^^ erg s LSI -i-6r303, located in die Northern hemi- 
sphere and inaccessible to HESS, has yet to be detected at TeV 
energies. While PSR B 1259-63 shares the property of being a 
TeV emitter with LS 5039, LSI -^61°303 and PSR B1259-63 
both have similar periodic radio outbursts that may be linked 
to both of them having a Be type companion. PSR B 1259-63 
has not been detected at GeV energies by EGRET (including 
at periastron passage, iTavani et all 1 199^ . perhaps because its 
long orbit only takes it close to the star for periods of time that 
are too limited for the sensitivity of the instrument. 

At closest, PSR B1259-63 is about 0.7 AU away from 
its companion, equivalent to the apastron orbital separation 
in LSI -i-61°303. The periastron separation of 0.1 AU in 



Table 1. Orbital parameters adopted for the y-ray binaries. 



System 


^"0* 


e 


M2 


ds 








(d) 




(Mo) 


(AU) 


(Re) 


(10" K) 


LS 5039 


3.9 


0.35 


23 


0.1-0.2 


9.3 


3.9 


LSI +6r303 


26.5 


0.72 


12 


0.1-0.7 


10 


2.2 


PSRB 1259-63 


1237 


0.87 


10 


0.7-10 


10 


2.7 



LSI -1-61° 303 is itself comparable to the orbit in the more com- 
pact LS 5039 (Table 0. There is a scaling between the three 
systems that, should they be governed by the same underlying 
processes, might lead to interesting comparisons. Indeed, their 
spectral energy distributions are similar in shape and luminos- 
ity (Fig.0, suggesting the same physics might be at work in 
these three known examples of "y-ray binaries". 

This work purports to investigate the pulsar wind scenario 
in light of the latest observations of these binaries. Reviewing 
current evidence, it is argued in §2 that there is as yet no rea- 
son to discard this possibihty in LS 5039 and LSI H-6r303. 
The scenario is then developed in the following sections. The 
basic assumptions are described in §3, showing how the phys- 
ical quantities scale with the parameters of the model. Simple 
expectations for the spectral energy distribution are derived in 
§4. A numerical model of the emission is elaborated in §5 and 
applied to the y-ray binaries in §6. It is shown that the radio 
emission from the shocked pulsar wind material can look like 
that of a compact jet emanating from an accreting microquasar. 
The conclusion set out in §7 is that the pulsar wind scenario 
provides a plausible, coherent and persuasive framework to in- 
terpret the emission from these y-ray binaries. 

2. Accretion or rotation-powered ? 

Several observations could securely tilt the balance towards 
the 'black hole jet', accretion-powered scenario or the 'pulsar 
wind', rotation-powered scenario in LS 5039 and LSI -i-61°303 
without detailed modelling. These are examined here in turn. 

2.1. Is the compact object a black hole ? 

A compact object mass greater than 3 M© would rule out 
the pulsar wind scenario. Constraints have been set in both 
LS 5039 and LSI -H6r303 by measuring the radial veloc- 
ity of the stellar companion. This enables the measure of the 
period, mass function, eccentricity and periastron angle. The 
mass function provides a strict lower limit to the mass of the 
compact object. In high-mass X-ray binaries the mass function 
is too small to provide a useful constraint: f(M ) = 0.053 Mp 
in LS 5039, f(M) = 0.011 Mg in LSI -h61°303 dCasares et alJ 
.2005 a b who have the best available datasets). 

A constraint on the compact object mass then requires 
additional knowledge on the mass of the stellar companion 
and the orbit inclination. Atmosphere fitting of the spectrum 
leads to an estimate for the mass (and radius) of the mas- 
sive star. Accuracy requires detailed stellar models and is in- 
herently difficult due to the high mass loss rates of early- 
type stars. Lines can originate at difi'erent depths and veloci- 
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ties in the strong (variable) stellar wind, leading to confusing 
measurements of radial velocities and difficultie s in spectral 
modelling (see e.g. discussion in 'B ahcail 11978V Parameters 
thereby derived have higher, probably irreducible, system- 
atics than t he quoted (statistical) errors might suggest (e.g. 
Table 2 of ICasares et all l2005b l). In fact, the orbital period 
of LSI +61°303 (26.5 days) cannot be determined from op- 
tical spectroscopy alone with high confidence and radial ve- 
locity studies use as an a priori the well-determined radio 
outburs t period (Hutchinas & Cramp ton, 198 1; Gregory. 2002: 
ICasares et alll2005ah . 

The parameter space derived bv lCasares et alJ (l2005a'b^ us- 
ing these techniques allows for both a neutron star or a black 
hole in LS 5039 and LSI -h61°303. The radial velocity am- 
plitude implies a low orbit inclination if the compact object 
is massive (/ * 30° for M > 3Mo in both LS 5039 and 
LSI -i-61°303), a good omen for a model with a relativistic jet 
pointed close to the line-of-sight. On the other hand, a 1 .4 
neutron star implies in both objects / ^ 60°. In LS 5039, 
the lack of X-ray eclipses (assuming point-like emission) by 
a companion of » 9 R© requires / < 65°, consistent with the 
above. Line broadening of the lines gives an upper limit on 
their rotation speeds which, since it must be lower than breakup 
speed, can be used to se t a lower lim it on the inclination for 
aligned spin/orbital axis. ICasares et al. ( 2005a,b) find i > 10° 
in both systems, again consistent with both a black hole or neu- 
tron star. Since the value of sin / is more likely to be high than 
low for a random distribution of inclinations, a neutron star 
would seem to be favoured over a black hole. 

The binary orbit of LS 5039 is quite compact for a high- 
mass X-ray binary (Porb ~ 3.9 days) and is only mildly ec- 
centric (e ~ 0.35). By comparison, LSI H-61°303 has Porb ~ 
26.5 days and e « 0.7. Orbital circularization and synchro- 
nization might be relatively ad vanced in LS 5039. This has 
prompted Casares et al.l (l2005bl) to use their measurement of 
rotational broadening to constrain the orbit in LS 5039 by as- 
suming that that the massive star is corotating at the angular 
orbital velocity at periastron. The compact object mass then 
becomes 2.7-5.0 M© with / = 24.9 + 2.8°, probably optimistic 
because of the possible systematics, but in practice ruling out a 
neutron star. 

This conclusion is not robust. Observations of the stellar 
companions in other high-mass X-ray binaries with similar 
periods do not systematically show corotation (Conti, 1978). 
Although synchronization is undoubtedly fast by astronomi- 
cal standards (some detached early-type binaries are observed 
to be circularized), short-lived wind-fed HMXBs do not have 
ages much greater than their typical synchronization timescale. 
These timescales are difficult to calculate precisely: in massive 
stars, they de pend on excited modes in the radiative envelope 
ll7ahn[ll97i . The age of the system can only be guessed at. 
Rotation-powered emission would probably imply an energetic 
pulsar with a short spindown timescale. Synchronization would 
then not have happened. 

In both LSI +6r303 and LS 5039, optical spectroscopy 
alone is unable to decide between a neutron star and black hole 
compact object. This work will assume a 1.4 M© neutron star. 



2.2. Is there accretion? 

Any sign of accretion on the compact object would dismiss 
the pulsar wind scenario. Accretion would have to be wind- 
fed: Roche lobe overflow as occurs in Cyg X-1 and some 
other HMXBs would entail much larger X-ray luminosities. 
The orbits and sizes of the massiv e stars in LS 5039 an d 
LSI -i-61°303 support this conclusion (Casar es et al l l2005albh . 
The rate at which mass is accreted from the stellar wind can 
be estimated. The Bondi capture radius is r^ = IGMdv^,, 
where v„, is the wind speed (about 2000 km s ') and Mc is 
the mass of the compact object. The Bondi mass accretion rate 
is then M x M„/(2ra/ds)^ where M„, is the stellar wind mass 
loss rate (about lO^^M© yr"') and c/, is the orbital separation. 
Using the LS 5039 orbital and stellar parameters in Table 1 
gives between 0.6-20 10^^ g s"' along the orbit in LS 5039, 
with a time-average of 5 lO'^ g s"'. With 10% radiative effi- 
ciency, this is too close to the sustained GeV emission at a level 
~ 10"^^ erg s In addition, although the X-ray luminosity is 
comparable to those of wind-fed X-ray binaries, the lightcurve 
does not show evidence for the large expected associated or - 
bital variability dRei g et all Eooi iBosch-Ramon et all l2005h . 
The slow, dense Be equatorial wind in LSI H-61°303 allows 
for higher accretion rates, in better agreement with the to- 
tal luminosity and the factor ~ 10 varia tions observed in the 
0.5-2 keV band lightcurve (less above. 'Ha rrison et alJl200flt 
iiartf & Paredes 1995, P aredes et al. 1997) . 

There are no signs of the large excursions in luminosity 
or spectra typical of (hard or soft) X-ray transients. The X-ray 
spectra are reminiscent of the low/hard X-ray state of LMXBs. 
There are some r adio/X-ray fluct u ations on timescales ~ hour 
in both objects llHarrisonetall l2000t iBosch-Ramon et all 
.2005) : but strong, fast variability on timescales of ~ second, 
expected from a lo w/hard state LMXB or accreting HMXB 
(Ivan de r 'Kiis','2006"), has not been detected ('Rib o et all [19991 
iHarrison et al._ 2000. although the low flux of the sources 
should be noted). The Thomson opacity is low in the stellar 
wind, hence variability might be smoothed out by scattering in 
a dense corona. There is no sign of such a corona in X-ray 
spectra. Indeed, the apparent lack of spectral cutoffs around 
either 30 keV or 100 keV (see Fig. 1) puts them apart from 
both accreting HMXBs and low state LMXBs. RXTE spec- 

tra of LS 5039 did show a strong, broad Fe Une dRibo et all 

I ' 

1999), but that could not be confirmed by other X-ray missions 

(Martpcchiaet al., 2005). LS 5039 Hes in the Galactic Plane 

and the Fe line is very Ukely due to Rid ge emission sampled 

by the large RXTE non-imaging aperture dBosch-Ramon et all 

am)- 

The radio and hard X-ray emission would be compatible 
with most of the emission coming fror n a non-thermal dis- 
tribution of particles in a rela tivistic jet (iRomero et al.L l2005t 
IParedes etalll2006t IPermer & Bottcheil l2006l) . Gamma-ravs 
would then be up-scattered photons. The puzzle is how such a 
jet could be fed by accretion without any of the usual signa- 
tures, most notably variability. In the y-ray binaries, the broad- 
band SEDs are stable even on timescales of years. 

The simplest hypothesis is that, as in PSR B 1259-63, ac- 
cretion does not occur in LS 5039 and LSI -i-6r303. If the 
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infall of material is stopped by pressure from a relativistic pul- 
sar wind, this will produce very little emission and variability. 
Instead, the emission will be due to the shocked pulsar wind. 
As in the relativistic jet scenario, a non-thermal distribution of 
particles provides the necessary ingredient to explain the SED. 
The spindown of the pulsar provides a stable energy source on 
human timescales. The termination shock of the contained pul- 
sar wind provides the conditions for particle acceleration. Only 
very mild variations are expected, due to fluctuations of the en- 
vironment (stellar wind), consistent with observations. 

2.3. Is the radio emission a microquasar jet? 

The detection of relativistic motion in the radio emission, like 
that of the microquasar GRS 1915+105, would secure the jet 
model. The radio morphology of LS 5039 and LSI -H6r303 
looks more like the compact jets seen in the low/hard X-ray 
state of low-mass X-ray binaries or Cyg X-1, with the signifi- 
cant difference that their radio spectral indices are steeper than 
the usual flat indices measured in the self-absorbed jets. The 
indices (Fy ~ v"*^^) rather hint at optically thin synchrotron 
emission from a dN oc E^^dE power-law distribution of elec- 
trons. 

The period of the radio outbursts of LSI -h61°3 03 is very 
stable and is used as a clock for the orbital motion ([Gregory, 
l2002h . The radio lumin osity of LS 5039 h as been steady at 
the same level for years JMarti et allll998ai) . The radio bright- 
ness asymmetry of LS 5039 and LSI -i-61°303 has been used to 
place constraints on (de)boosting of the putative relativistic jet . 
This yields a mild bulk fl ow velocity of 0.2-0. 3c jParedes et all 
l200rt'Massi et alll2004 ). There is no evidence in either system 
for highly relativistic bulk motion directed more or less towards 
the line-of-sight; the main argument favouring a jet interpreta- 
tion is based on the morphology of the radio emission. 

As mentioned in §1, fast-moving young and old millisec- 
ond pulsars leave trails of X-ray, optical and radio emission 
that can reach parsec lengths, and whose properties scale with 
the distance where the wind terminates due to conta i nment 
by the surrounding in terstellar medium ( I Wang et all Il993t 
lArons & Tava ni'. '1993"). Even the collision of standard stellar 
winds can pr oduce detectable large- scale, jet-looking, radio 
emission (e.g. IPoughertv et"all 120031 . The radio emission of 
LS 5039 and LSI -h61°303, seems to reach scales of up to sev- 
eral 100 AUs ( Paredesetal.L 2002: Massi et al.. 2004.) , double- 
sided in the case of LS 5039. Energetically speaking, this emis- 
sion can easily be a pulsar nebula. However, isolated ms pulsars 
have time to develop long (one-sided) nebula away from their 
direction of proper motion. In the y-ray binaries, the direction 
of such a nebula would be affected on short timescales by or- 
bital motion. Radio observations of LSI +61°303 have shown 
curved radio emission on a scale of 50 mas, interpreted as due 
to precession of a jet as in SS 433 ('Massi et all 1200 4*). Could 
this morphology instead be a result of the pulsar orbital mo- 
tion? 

In this work, it will be argued that the resolved radio emis- 
sion in LS 5039 and LSI +61°303 is not due to a relativistic 
jet akin to those of AGN or microquasars (and implying accre- 



tion), but arises from shocked pulsar wind material outflowing 
from the binary. 

2.4. Can the compact object be a pulsar? 

The X-ray and radio properties of LS 5039 and LSI -H6r303 
sets them apart from other X-ray binaries, with no conclusive 
signs of accretion going on. At the same time, the luminosity, 
hard X-ray power law extending to high energies, weak vari- 
ability and optically thin radio emission are typical for emis- 
sion from young pulsar s with a high spindown power (e.g. 
iGaensler & Slan e. 2006). Examining the pulsar wind hypoth- 
esis therefore seems justified. 

In order to avoid disc or magnetospheric accretion, the rel- 
ativistic wind of the pulsar must be able to quench the in fall of 
stellar matter attracted to it JlUarionov & Sunyaevlll97^ . If all 
of the pulsar spindown power E is transferred to this relativistic 
wind, then writing the pressures gives E > AM\\^.c where M is 
the Bondi rate calculated above. Writing the spindown power 
as a function of the pulsar magnetic field B and period P, then 
P < 230B[f ms, with M15 = MjlQ^^ g s"' . A young, ms 
pulsar would push out any settling accretion flow. 

The young age is not problematic. Radio searches have not 
concl usively identified an association with a supernova rem - 
nant Jprail et all fl987t iMartf et all Il998at iRTbo et all l2002l) : 
but neither have they been able to do so for PSR B1259- 
63, which has a measured spin down age of 3 10^ years 
Jjohnston et alIll992h . lRib6 et alJ ( 12002 ) put an upper limit of 
about a Myr to the age of LS 5039 by tracing back the proper 
motion of the system to the plane. Nitrogen enrichment in the 
companion's atmosphere also suggests a young age in LS 5039 
(McSwain et al. 2004; rotational-mixing could also explain the 
enrichment. lCasares et alJ2005bl) . 

Detecting the radio pulse would nail down this scenario. 
Unfortunately, free-free absorption in the dense stellar winds 
of LS 5039 and LSI -H6r303 would suppress flie signal. 
Using the Rosseland mean opacity, the radial optical depth 
at 1 GHz from a point at a distance d from the star is r 
3 ■ ^^^Kj^timJA'^^^cnAA for a 10"^ yr"' coasting 
wind with a speed 2000 km s"' and temperature 10^ K; is 
taken to be 0.1 AU, about the periastron orbital separation in 
LS 5039 and LSI -H6r303. High opac ities are also found in 
LSI -H6r303 (iTavlor & Gregorylll982h . In PSR B1259-63, the 
suppression of pulsed emission is observed close to periastron 
Jjohnston et alIll992tlMelatos et allll995l) . 

Pulsed emission in the X-ray/GeV band would serve just 
as well. Such a detection does not seem realistic without prior 
knowledge on the pulse period derivatives given the sensitiv- 
ities achievable even by GLAST. In addition, the pulsed GeV 
signal in LS 5039 may be smoothed out by emission from the 
e^e^ cas cade initiated by absorption of TeV y-rays on stellar 
photons JPubuslEoO ^. 

As the pulsar rotation gradually decreases, the y-ray bi- 
naries would become standard accretion-powered high mass 
X-ray binaries, perhaps brighter in the X-ray sky but proba- 
bly without significant emission at y-ray energies. That three 
examples of such short-lived phases may have been found 
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perhaps implies a large population of descendant HMXBs 
of which some might be the enshrouded hard X-ray sources 
discovered by INTEGRAL. Some of the unidentified EGRET 
sources located in the Galactic Plane may also turn out to be 
pulsars in binar ies. Massive X-ray binary popul ation synthesis 
calculations bv lMeurs & van den Heuvell ill 989 ) predict about 
30 binaries in the brief young pulsar stage in the Galaxy, con- 
sistent with three examples being found within 3 kpc. 

The conclusion is that LS 5039 and LSI -H6r303 could 
very well harbour a young pulsar with a strong relativistic wind 
like that of PSR B 1259-63. A model based on this scenario is 
explored in the next sections. 



3. Scaling relationships in the binary plerion 

The scenario, identical to that proposed for PSR B 1259-63, 
is the following: a young pulsar with a spindown rate of or- 
der 10^^ erg s"^ g enerates a strong rel ativistic wind beyond 
the hght cylinder jRees & GunnL [197I . The pulsar wind is 
assumed to be radial, isotropic and composed of principally 
mono-energetic electrons/positrons with a Lorentz factor jp. 
All the spindown power is transferred to the leptons, with a 
small fraction going to the magnetic field. The ratio of mag- 
netic to kinetic energy cr is a parameter of the model. 

The pulsar wind is contained by the stellar wind (in a 
plerion like the Crab, the container is the supernova remnant 
material). A coUisionless shock forms beyond which the lep- 
tons from the pulsar wind are accelerated and isotropized. The 
shocked pulsar wind material is separated from the stellar wind 
by a contact discontinuity across which mixing might occur 
(neglected here). The shock has a 'bow' or 'comet' shape with 
a tail extending away from the stellar companion. The radio 
emission occurs in this nebula. Studying such a complex inter- 
action is simplified by noting everything scales with the stand- 
off distance (see e.g. Bucciantini et al., 2005). The standoff dis- 
tance Rs is the radial distance from the pulsar to the star at 
which the ram pressures equilibrate (§3.2). 

Most emission occurs in the vicinity of Rs, whose value can 
be calculated given the orbital parameters, and the stellar wind 
density and velocity. Straightforward estimates indicate that the 
physical conditions at the shock are set by only three quantities: 
the shock distance R^ from the pulsar, the orbital separation ds 
and the combination crE, which describes the pulsar wind en- 
ergy and magnetization. The overall properties of the resulting 
spectral energy distribution are then deduced in the next sec- 
tion (§4). A quantitative description including emission from 
the nebula is given in §5-6. 



3.1. The pulsar wind 

The pulsar wind is described by the Lorentz factor jp of the 
electrons/positrons, the spindown power E and the ratio of 
magnetic to kinetic energy cr. The (proper frame) density iii 
and (observer frame) magnetic field Bi in the pulsar wind up- 
stream of the shock are then derived from (.Kennel & CoronitL 
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Ann\y^pmeC^ 
1 -1-0- 

4;r 



(1) 
(2) 



10^^ erg s ' and the pulsar mag- 



Most of the energy could actually be carried by ions, which 
are probably a necessary ingredient to accelerate the leptons 
at the termination shock ( Arons & Tavani, 1994). Here, the ki- 
netic energy of the wind is assumed to be only in the elec- 
trons/positrons. 

In PSR B 1259-63 £ is 8 
netic field B is about 3 ■ 10" G (corresponding to a spin period 
P — Al .1 ms and period derivative P - 2.3 ■ 10"'^). Thereafter, 
E wiU be expected to be of order 10"^* erg s Models of 
the Crab plerion yie ld values of order 10~^ for cr and y„ in 
the 10 ^ - 10^ range jRees & Gunnl Il974t iKennel & CoronitH 
Il984ah . Explaining the implied conversion of electromagnetic 
energy (close to the pulsar) to kinetic energy (far from the pul- 
sar) is a key question of pulsar physics. Values for cr could be 
higher if the shock occurs clos er to the pulsar than in C rab, as 
proposed for PSR B1259-63 bv lTavani & Aronsid 19971) : in any 
case cr will be expected smaller than 1. 

3.2. The standoff distance 

The standoff distance point R^ of the pulsar wind is given by 



AttR^a 



(3) 



where Vy, is the pulsar radial orbital speed calculated from the 
binary parameters, p^. is the stellar wind density and Vv, the 
stellar wind speed at the standoff distance. An estimate of Rs 
can be obtained by balancing the pressures of the two flows, 
assuming a stellar wind coasting at a speed Voo and a negligible 
orbital speed: 



Rs 



ds 



1 + (Pn./E) 



1/2 



(4) 



where P„,/c = MwVco is the stellar wind momentum flux with 
M„, the stellar mass loss rate, supposed isotropic here, and ds is 
the orbital separation. To fix ideas, the radiatively-driven wind 
of LS 5039 has M„, ^ 10"^ Mq yr^' and v^, ^ 2000 km s"' 
which gives P„ ^ 5 ■ lO" erg s"' ilMcSwain et al.[Eo04l) . When 
the stellar wind is more powerful than the pulsar wind, in the 
sense that P„ » E, then the standoff point can be written as 



Rs 



1.5 



lQ''cmdo.iP;^f,E\l' 



(5) 



where do,i = ds/O.l AU is the orbital separation (0.1 AU at pe- 
riastron in LS 5039 and LSI -H6r303), £36 = E/IO^^ erg s"' 
and Ph.38 = M„,Vcoc/10-'^ erg s"'. The standoff distance stays 
at about a tenth of the binary separation throughout the orbit. 
Writing down explicitly the dependence on the stellar wind 
(Eq. 13 is easy using Voo and P,,,. Eq. 0]is used for the slow 
equatorial winds of Be stars (§6). In the following, it proved 
more convenient to explicitely keep Rs as a parameter, derived 
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under the plausible assumptions described above. Some of the 
uncertainties in the geometry of the shock region or the stellar 
wind parameters (particularly in the wind acceleration region) 
are thereby hidden in Rs. 

3.3. Downstream flow at standoff 

Downstream values for the magnetic field and particle density 
are calculated using perpendicular MHD shock jump condi- 
tions at Rs (Kennel & Coro nitL ■1984a.) . The post-shock mag- 
netic field is (for cr <s; 1) 



B = 3(1 -4o-) 



E/c cr 
/?? l+o- 



1/2 



5 (E^ecTiy'^Ril G 



(6) 



where 0-3 = cr/ 10"^ and/; 11 =«,/10"cm. 

The downstream particles are assumed to be accelerated to 
a canonical dny oc y^^dy power law from a minimum energy 
ymin to a maximum energy ymax- For a given y^ax, the normal- 
isation and lower limit ymin of the post-shock particle distri- 
bution are found by matching to the downstream particle and 
ene rgy densities. For a 7'^ pow er law and cr <s; 1, is given 
by jKennel & Coronitiill984bl) 



,ln 



V2 



(7) 



For a post-shock distribution covering 2-3 decades in energy, 

Tmin ~ O.lyp. 

3.4. Timescales at Standoff 

Particles are assumed to be accelerated efficiently up to the 
maximum energy y^^^ above which the timescale at which par- 
ticles radiate their energy becomes shorter than the acceleration 
timescale. The gyroradius is used to give an estimate of the ac- 
celeration timescale to a Lorentz factor 7 



ym^c 
eB 



0.01 76 (£'360"3) ^^^Ru seconds 



(8) 



using the downstream magnetic field (Eq.|6J and writing 75 = 
7/10^. For ion dominated shocks (.Arons & TavanL.1994.) , the 
relevant timescale is given by the ion gyroradius and the max- 
imum electron energy is ~ ypinp/nig (jp is then the upstream 
ion Lorentz factor). The maximum 7 reached must have a gy- 
roradius smaller than the characteristic scale of the shock R^, 
implying 7 < 3 ■ 10^ (£'36cr3)'^^. Particles do not reach such 
high 7 as they then lose energy to radiation more rapidly than 
they can be accelerated. 

The synchrotron timescale is 



9mlc^ 1 
4 e^y B^ 



30 7g' (£360-3) '/Jjj seconds 



so the maximum 7 such that fj 



sync 



> he 



5 10' (£360-3) 



-1/4 r>l/2 
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(9) 



(10) 



which, in principle, could result in a maximum photon en- 
ergy of 25 TeV. Photons up to energie s > 4 TeV are de- 
tected by HESS (fAharonian et all Eo05al) . requiring electrons 



with Lorentz factors high enough that 7max'MeC^ ^ 4 TeV, i.e. 
Tmax ~ 10^. This is below the limit from synchrotron losses 
(Eq.llO> for all plausible values of the pulsar o--parameter. 

The inverse Comp ton timescale in the Kle in-Nishina 
regime yhv » mgC^ is jBlumenthal & Gouldlll97nl) 



fKN -5.2 76 



3eV 
'kfZ 



In 



4ykT^ 



1.98 



seconds(l 1) 



where the stellar (blackbody) radiation density has been low- 
ered by geometric dilution and Rs is assumed to be small com- 
pared to the orbital separation ds- The stellar temperature 
varies between 2^-10'* K with the stellar radius » IQRq for 
LS 5039, LSI +61°303 and PSR B1259-63 so that 



fKN =20 764 Jin 76-1- 1-3]- 



7^* 4^*^ioseconds 



(12) 



with r*,4 = r*/40,000 K and R^^io 
Thomson limit yhv <K m^c^ then 



Ri^/lORo. In the 



9 mjc^ 
tj - - 



1 



4 SttJ/* 



= 40 73'4i7'*,4«* 



10 



seconds 



(13) 



where f/* = {o-BT*/c){Ri,/ds)^ is the stellar photon density. 
Most of the stellar photons have an energy of UkT-i, ^ 9 eV 
so that the turning point between Thomson and Klein-Nishina 
regime is for 7kn ~ 511 keV/9 eVa; 6 ■ lO'*. The interaction 
timescale is then at a minimum of a 7 s (see Fig.^J- 

The timescale for particles in the pulsar wind to reach the 
standoff point is 



Rs/c ^ 3 Rii seconds 



(14) 



which is only slightly shorter than the minimum inverse 
Compton timescale derived above. This is important since if the 
pulsar wind is composed of mono-energetic particles with yp a 
6 ■ lO'*, then a fraction of the spindown power can be radiated 
by bulk inverse Compton interactions in 7-rays before reaching 
the termination point (Compton drag). There may then not even 
be a shock. (The upstream magnetic field is frozen in the MHD 
flow so upstream particles do not radiate synchrotron.) The re- 
sulting 7-ray spectrum would displ ay line-like emission with 
a strong orbital phase dependence ("B ogovalov & Aharonianl 
1^000; Ball & Kirk, 2000). This possibility appears excluded 
by the power-law HESS spectra extending to a few TeV 
in LS 5039 and PSR B1259-63 (lAharonian et all l2005albh . 
Compton losses are weaker in LSI -H6r303 and PSR B 1259-63 
(smaller and larger separations) than in LS 5039. 

Downstream of the shock, the inverse Compton timescale 
is greater than the synchrotron timescale when 7 < 7s where 



7s * 1.2- lO*- (Rn/dQ.i)(E,,cr3r"^T^,4R 



*,io 



(15) 



using Eq.OandEl The maximum electron energy ymax is in- 
deed set by synchrotron losses (Eq. II Oi l, and not by inverse 
Compton losses. Perhaps counter-intuitively, high energy par- 
ticles lose their energy mostly through synchrotron radiation, 
although the radiation density at Rs from the luminous stel- 
lar companion is much larger than the magnetic field energy 
density (U^/Ug ~ lO\Rn/doA)\E3(,o-i)-^). The difference 
in densities is compensated by the reduction of the inverse 
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Compton cross-section in the Klein-Nishina regime. Note that 
js is independent of the orbital separation when the stellar wind 
dominates (R^ oc d^, see Eq.|5}- 

When cr <K 1, the post-shock material is found from 
the jump conditions to i nitially flow away at v x c/3 
jKennel & Coronitilll984al) so the escape timescale is of order 

tesc ~ — rr ~ 10 7^11 seconds (16) 

c/3 

Here, '/aia ~ O.ly^ > 7kn so particles radiate in the Klein- 
Nishina regime up to ys- Their radiative timescale tic ~ 7- 
20 s is comparable to their escape timescale. Particles with y > 
7s lose energy to synchrotron radiation, on a timescale fsync $ 
30 s that is, again, comparable to their escape timescale. The 
synchrotron radiative efliciency increases with decreasing Rs 

since Tsync/^esc ^ Rs- 

4. Expected high energy spectrum 

The timescales given above are sufficient to get an approxi- 
mate idea of the expected high energy emission from the pulsar 
wind shock region. The aim is to provide a simple framework to 
understand the results from the more complex modelling pre- 
sented in §5-6. 

4.1. Spectral shape 

Since the radiation timescales are comparable to or (much) 
smaller than the flow timescale, the situation can be approxi- 
mated to steady-state injection of power-law electrons interact- 
ing in the Klein-Nishina regime with an isotropic bath of quasi 
mono-energetic photons. Mode rski et alJ ( l2005l) have studied 
exactly this case, taking into account synchrotron emission in a 
uniform magnetic field. The overall characteristics of the spec- 
trum emitted in the vicinity of the pulsar can be deduced from 
their results and the calculations above. If particle cooling is not 
fast enough, then particles will travel to regions of lower energy 
densities, modifying the spectral energy distribution from that 
presented here. The steady state assumption is reasonable for 
the high energy spectrum (§5). 

[Moderski et alJ ll2005l) show that the steady-state electron 
distribution displays three regimes: (a) at low energies, inter- 
actions in the Thomson limit steepen, as usual, the distribu- 
tion by 1 in index; (b) at high energies, synchrotron emission 
also steepens the distribution by 1; (c) at intermediate ener- 
gies, where inverse Compton interactions in the Klein-Nishina 
regime dominate, the decrease in cross-section is such that the 
distribution index is unchanged or hardened at most by < 0.5. 
The hardening decreases with the ratio U^/Ub- The change be- 
tween regimes (a) and (b) occurs at the transition to the Klein- 
Nishina regime yKN- Case (a) does not happen here since we 
initially have y^nin > Tkn- The change between regimes (b) 
and (c) occurs at ys where synchrotron losses start to domi- 
nate. A dN ~ y^^dy injection spectrum will therefore lead to a 
steady-state electron distribution going from dN ~ y^^-^dy to 
y^^ above ys- 
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Using the downstream magnetic field (Eq. |6j, the charac- 
teristic synchrotron frequency for an electron with an energy 

ym^c^ is 

/;y^ 100y2(£36(^3)'^Xi' keV (17) 

The three frequencies relevant to the spectral shape, corre- 
spo nding to ymin, ys and y^^x are (using y^in ~ O.ly^ and 
Eqs. FToFTsI : 

/iv.sy„c,m,„ ^ 1 y2 6(£36cr3)'^'/?„' keV (18) 
hVsyncS ^ l50(E36CT3)-"\T^,4R*,w)\Ru/4i) keV (19) 
/jVsyncmax ~ 250 Me V (20) 

The expected synchrotron spectrum jModerski et allEo05l) is a 
hard power law vFy ~ v"^ from about 1 keV initially (ymin) up 
to a break energy » 100 keV (ys). The spectrum then flattens, 
with vFy ~ constant up to s; 250 MeV (ymax)- Since R Jds is 
almost constant (Eq. 5), the break frequency changes as dj^ 
along the orbit because of the variation of the magnetic field 
at the standoff point. The maximum synchrotron frequency is 
independent of all parameters. 

Inverse Compton interactions occur in the Klein-Nishina 
regime so that the characteristic frequencies are given by hv * 
ynigC^ 

*50yp,6GeV (21) 
/2v,omp,s * 0.6(E,6(r3r"\T^AR*,io)(Ru/doA) TeV (22) 
hv,o„, ^ 25 (E36cr3r"^R\{^ TeV (23) 

iModerski et all tOO^ find the inverse Compton spectrum for 
a y"^ injection is constant or rising slightly in vFy up to Vcomp,s 
(their Fig. 6). Above an energy y^nieC^, the electron energy dis- 
tribution steepens to y"^. The spectral index steepens by ^ 1.5 
(their Eq. 31) and vFy ~ v"'-^ is expected above 0.6 TeV. 

Note that photons emitted with energies above about 
30 GeV can interact with stellar photons in the binary sys - 
tem to produce pairs yy e^e' llGould & SchredeS Il967h . 
This leads at least to a hardening (softening) of the observed 
VHE y-ray spectrum above (below) a few 100 GeV jPubusl 
'2006'). This can produce modifications to the electron distribu- 
tion since pairs are created (this is further discussed in §6.1.1). 

4.2. Luminosity 

Assuming the situation is steady-state and that the pulsar spin- 
down power is largely given to the particles (cr <s; 1), the syn- 
chrotron and inverse Compton luminosities should be, to order 
unity, equal to E. With a y"^ injection, the inve rse Compton and 
synchrotron luminosities are roughly equal jModerski et all 
.200S) . 

With such short radiative timescales, only a moderate spin- 
down power of lO-'^ erg s"' is sufficient to achieve a high ra- 
diative efficiency in X-rays and y-rays. This has little to do 
with the assumptions made on the pulsar wind y^, cr or £ (if 
the wind is purely composed of electrons and positrons: a pul- 
sar wind with a significant ion fraction would require a higher 
E to obtain the same luminosity). The efficiency comes from 
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the high radiation and magnetic energy densities at the stand- 
off point which allow fast cooling of all of the injected parti- 
cles. This contrasts with pulsars interacting with their super- 
nova remnant or the interstellar medium, whose termination 
shock occurs much farther out, at correspondingly lower en- 
ergy densities. 

The main factor setting the high energy luminosity is the 
shock distance Rs or, alternatively, the orbital separation ds 
(since R, oz ds in Eq.|5}- The shock distance enters as in the 
synchrotron timescale compared to Rs in the escape timescale. 
The smaller Rs is, the more synchrotron emission dominates 
and the more the pulsar can be expected to be radiatively effi- 
cient. However, this is at the expense of the inverse Compton 
flux since the associated timescale stays constant with Rs. 

The dominant radiation process does not change along the 
orbit for a given electron energy (because ys ~ Rs/ds), but the 
probability that this electron can escape from the region of high 
Ui, and Ub without significant radiation increases with ds (or 
Rs). The steady-state assumption then breaks down. Conditions 
change as particles flow away from the shock region. The over- 
all nebula must then be taken into account to model the emis- 
sion (§5). 



that have cooled and escaped from the vicinity of the pulsar 
(see §5-6). 

The GeV emission may be problematic as the model pre- 
dicts a gap between 250 MeV and 50 GeV when jp = 10^. 
There is no room for increasing the synchrotron emission so as 
to include at least a fraction of the EGRET band: the maximum 
synchrotron frequency is fixed by the assumptions on ymax- 
Advocating higher y^ax seems unreasonable. Moreover, for 
large orbital separations as in PSR B1259-63 or LSI -h61°303 
(at apastron), the maximum electron energy is limited by 
the Larmor radius (supposed < Rs) rather than by radiative 
losses. These range from 2 MeV (apastron, polar wind) to 80 
MeV (periastron, equatorial wind) for PSR B 1259-63 and from 
30 MeV (apastron, polar wind) to 250 MeV (periastron, equa- 
torial wind) for LSI -i-61°303. The easiest solution consists in 
lowering the Lorentz factor jp of the particles in the pulsar 
wind to 10^. This lowers Vcomp.min to a few GeV whilst leaving 
the spectral breaks unchanged. The inverse Compton timescale 
becomes long compared to fesc at the low energy end of the 
electron distribution and, as for radio, high energy emission 
between 250 MeV and a GeV will also come from escaping 
particles. 



4.3. Comparison to SEDs 

The characteristic frequencies for all three y-ray binaries were 
calculated using the parameters in Table ^ and additional as- 
sumptions on their stellar winds. For LS 5039, the wind is ra- 
diative with M„, = 10"^ Mq yr"' and Voo - 2400 km s"'. For 
LSI -h61°303 and PSR B1259-63 two winds were assumed: a 
dense, equatorial outflow and a fast, polar wind similar to that 
in LS 5039. Details are given in §6. The standoff distances Rs 
and frequencies are given in Table |2 

Anticipating on §6, the pulsar was assumed to have jp = 
10^, cr = 10"^ and E - 10^^ erg s"' in all three systems. 
Note that E and cr always appear in the same proportion so 
that for given Rs there are really only two parameters, jp and 
(tE. The characteristic frequencies are not very sensitive to crE. 
Interestingly, the SED is also indifferent to the exact value of 
jp but for the low energy end of the inverse Compton spectrum. 

The expected breaks match reasonably well the observed 
SEDs of LS 5039, LSI -H6r303 and PSR B1259-63, notably 
the rising then flattening X-ray to y-ray spectra. An important 
observable feature is that, all other parameters being equal, the 
break frequency in X-rays ~ l/ds should be lower at apas- 
tron than at periastron. At higher energies, the break Vcomp.s in 
LSI -h61°303 and PSR B1259-63 occurs at frequencies a fac- 
tor 2 lower than in LS 5039, because of the lower (TableQ). 
The TeV spectrum of PSR B 1259-63 is indeed steeper than that 
of LS 5039. Unlike what happens with synchrotron, the spec- 
tral break in the inverse Compton emission should not change 
along the orbit if Rs is proportional to ds (Eq. 5). 

However, there are two shortcomings: the radio emission 
and the EGRET detections of LS 5039 and LSI -i-61°303 at 
GeV energies (Fig.^. The radio emission cannot be addressed 
by the present scheme: the minimum synchrotron frequency is 
ss 1 keV and radiation well below must originate from particles 



5. A numerical model of the pulsar nebula 

The detailed spectral energy distribution, including evolution 
of the physical conditions in the flow, is calculated here with 
the aim of accounting for the resolved radio emission as well as 
providing more quantitative comparisons. The assumptions are 
described in this section, with LS 5039 used as an illustrative 
example when necessary. The application and comparison with 
observations are made in §6. 



5.1. Evolution of tlie particle distribution 

Particles are continuously injected with a y"^ power-law energy 
spectrum between ymin and ymax beyond the termination shock, 
in an emitting region of characteristic size Rs- Other distribu- 
tions have not been considered for the sake of simplicity. The 
downstream particle densit y, in the limit of small cr is given by 
dKennel & Coronitilll984bi) 



n2 ~ SniYp X 980 y^j '^ae/^if particles cm ^ 



(24) 



using Eq. |2 The normalisation of the injected particle power- 
law Qiy) is then found from 



f 



Q{y)dy * -njnR"; ^ 3 ■ 10^^ Te'^'sfe particles s"' (25) 



Particles move downstream away from the pulsar with a 
speed V (initially a; c/3) whilst emitting radiation. The resulting 
'cometary' nebula is modelled by dividing it into small (loga- 
rithmic) subsections where the physical parameters (magnetic 
field, density etc) are homogeneous. The nebula is assumed 
to expand conically. The evolution of the distribution function 
N{y, z) along the nebula is given by 



— - — 

dz By 



dz 



+ Q(y)S(z - 1) 



(26) 
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Table 2. Characteristic frequencies of the expected spectral energy distributions. 
The standoff distance (first number) is calculated at periastron and apastron us- 
ing Eq. I^and the stellar wind parameters in §6. Both a polar and an equatorial 
wind are considered for the Be stars. The break frequencies in the synchrotron and 
Compton spectra are then calculated using Eq.^]and Eq.|22]in ^4.11 The pulsar 
wind parameters are the same for all the objects: E - Iff'^ erg s"', cr - 0.01 and 
jp - 10^. The star temperatures are such that inverse Compton interactions are in 
the Klein-Nishina regime. 



LS 5039 




periastron = 0. 1 AU 




apastron d, = 0.2 AU 




2 


10" cm 100 keV 0.4 TeV 


4 


10" cm 200 keV 0.4 TeV 


LSI -(-61° 303 




periastron =0.1 AU 




apastron d^ = 0.7 AU 


polar 


8 


10" cm 100 keV 0.8 TeV 


4 


10'2 cm 15 keV 0.5 TeV 


equatorial 


5 


10'" cm 6 keV 0.05 TeV 


2 


10'2 cm 5 keV 0.3 TeV 


PSRB 1259-63 




periastron = 0.7 AU 




apastron d^ = 10 AU 


polar 


3 


10" cm 12keV 0.5 TeV 


4 


10" cm 1 keV 0.5 TeV 


equatorial 




10'^ cm 4 keV 0.2 TeV 




not applicable 



where z is the distance along the nebula normalised to the 
standoff distance R,. The radiative and adiabatic losses for an 
electron (positron) are: 



dy 
dj 



V 



'dy 
, df 



sync 



(27) 



The synchrotron radiation loss is the usual ysync ~ y^B^. The 
ene rgv lo s s to in verse Compton radiation is calculated using 
the Llonesl ill 9681) kernel, which is a continuous approximation 
although electrons can lose significant amounts of approxima- 
tion in discrete encounters with photons in the Klein-Nishina 
regime. Furthermore, the Jones kernel supposes an isotropic 
photon distribution. This should not be too bad an approxima- 
tion as the angle-dependence of the cross-section is limited in 
the Klein-Nishina regime. 

The adiabatic losses for two dimensional expansion of rel- 
ativistic material are 



dlnn'^^^ 
dz 



2y V 



(28) 



The evolution of the physical parameters of the flow (B, 
V etc) with z is described in the next section. Once these are 
determined, the differential equation governing the evolution 
(Eq. I26> is numerically integrated along z u sing a logarithmic 
grid in y, following the scheme proposed bv lChang & CoopeJ 
The results presented thereafter typically use at least 
1000 points for the particle energy grid and 200 points in fre- 
quency V to represent emission spectra. The results have been 
checked to be independent of the numerical details. 

5.2. Evolution of physical conditions in the flow 

The evolution of the nebular flow speed v and magnetic field B 
with z are computed by using the relativistic Bernouilli equa- 
tion stating that 



B^ 



-I- 



+ P + ne \ - constant 



along streamlines for stationary, adiabatic motion. F is the bulk 
Lorentz factor of the flow F = (1 - v^/c^)"'^^. The addi- 
tional relationships for the magnetic field B, density «, pres- 
sure P and internal energy e come from the conservation of 
(toroidal) magnetic flux Bvz/F, of the number of particles nvz^ 
and the relativistic equation of state P - n{e - mc^)/3 oc n"*/^ 
dKennel & Coronitilll984ai) . 

Close to the shock, the total energy is dominated by pres- 
sure. The density is then roughly constant, the flow speed de- 
creases as V oc l/z and the magnetic field increases because 
of flux conservation (Fig. 3). The flow becomes magnetically 
dominated at large radii (z > 10) at which point the density 
varies as 1/z^, the magnetic field decreases as B oc l/z and the 
asymptotic flow speed is crc (Kennel & Coroniti, 1984a). 

The present approach is admittedly simplistic. MHD simu- 
lations of a spherically symmetric pulsar wi nd interacting with 
a constant external medium exist (Buccianti ni et alll2005l) and 
give insights into possible shortcomings. The shape of the pul- 
sar termination shock is bullet-shaped rather than spherical. An 
asymmetric pulsar wind would surely change the aspect of the 
flow (as it does with the Crab and Vela pulsars). Shocked ma- 
terial flows along two channels, a fast one with material from 
the forward termination shock and a slower one with material 
from the backward facing termination shock. The nebular flow 
loses energy to radiation, contradicting the adiabatic approx- 
imation. Collimation shocks may re-energise particles along 
the jet; mixing can occur at the surface discontinuity between 
the fast-moving shocked pulsar wind and the slow-moving stel- 
lar wind. Mass loading will slow down if not disrupt the flow. 
At larger distances, the flow shears because of orbital motion 
(§5.4), leading to additional complications. Finally, the flow 
merges with the stellar wind and/or ISM on scales of ~ 0.01 pc 
based on the densities. 

5.3. Calculated spectrum 



Neglecting any effect of the orbital motion (see ^5.4> . the over- 
all spectrum results from the summed contribution of the syn- 
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Fig. 3. From left to right: evolution in the pulsar nebula of the flow speed, magnetic field and density as a function of the 
normalised distance z- The rightmost panel shows the time needed to reach a given distance z in the flow. At z » 800 id - 
zRs ~ 10 AU), the elapsed time is equal to the orbital period (dotted line). The model parameters are E = 10^^ erg s y = 10^, 
cr - 0.01 and Rs - 2 ■ 10" cm, appropriate for LS 5039 at periastron. 



chrotron and inverse Compton spectra (on stellar photons) cal- 
culated for subsections dz of the nebula extending away from 
the star. As was done for energy losses, the particle distribu- 
tion and radiation fields are assumed spatially homogeneous 
and isotropic in each subsection. The inverse Compton spec- 
trum takes into account the reduction of the cross-section for 
interactions in the Klein-Nishina regime via the Jones kernel. 
The radiation field of the stellar companion is simply diluted 
by the inverse square of the distance to the star (11* cc d^^ with 
d* = d, + (z- l)Rs). 

Fig -Ejillustrates the evolution of the emission spectrum and 
particle distribution with increasing distance z to the pulsar in 
the case of LS 5039. The first 3-4 spectra show the initial evolu- 
tion towards the steady-state described in §4. The particle dis- 
tribution above ys ~ 10^ cools rapidly at high energies due 
to synchrotron losses, resulting in a plateau down to vs as ex- 
pected. Evolution then sets in around z ~ 10, with the spectral 
luminosity decreasing with B. The particle distribution steep- 
ens with emission sweeping the optical to IR bands. Further 
evolution, corresponding to peak emission in mm to radio, is 
set by adiabatic losses. 

Self-Compton radiation is negligible. Synchrotron emis- 
sion was optically thin down to low radio frequencies in all 
the models computed for §6. This is a consequence of the flow 
geometries, which are discussed in ^5.41 Still, synchrotron self- 
absorption might be possible in a flow directed predominantly 
along the line-of-sight. 

In the absence of orbital motion, the spectrum seen by the 
observer is simply the integrated emission along the 'comet' 
tail nebula, just as for isolated pulsar interacting with the ISM. 
However, as can be deduced from Fig. 0] the radio emission 
originates from regions several AU away from the pulsar (z > 
100). The time At it takes for particles to reach this location, 
about a day here, is a significant fraction of f orb = 3.9 days. 
Orbital motion will have changed the inner boundary condi- 
tion so the spectrum seen by an observer will be a composite 
of SEDs from all orbital phases 0. Here, the only parameter 
dependent on the orbital phase is the pulsar standoff point 
Other dependences can be envisioned if, for instance, the pulsar 
wind is anisotropic or if the magnetization parameter is a func- 



tion of distance to the light cylinder (with higher cr expected 
close to the pulsar). 

The relevant timescales to compare are the flow timescale 
^flow = d/vto reach a distance d from the shock and the orbital 
motion timescale Torb - t/j/vorb- The orbital timescale can vary 
greatly for a highly eccentric orbit. Nevertheless, taking v equal 
to its asymptotic value crc, the distance where orbital motion 
becomes important will on average be 

dn ~ o-cPoib/27r (30) 

and the associated timescale, Poi-b/ln. As it turns out, the syn- 
chrotron emission from LS 5039 depends little on the orbital 
phase because R Jds is constant with (p. The spectrum seen by 
an observer should therefore be close to the integrated spectrum 
shown in Fig. Eleven at low frequencies (see ^6.1> . 

5.4. Appearance of the cometary nebula 

On a small scale, the shocked material flows away from the 
binary on a straight path, following the direction given by the 
vector difference of the stellar wind and orbital speeds Vw-Vorb- 
On a larger scale, the material 'bends' as the orbital motion 
becomes important at ~ £/„ (about 1 AU for LS 5039; Ea. l30> . 
Parcels of material are still on outward straight paths but the 
orbital motion shears the nebula. The opening angle of the flow 
(~ Rs/ds) is small so that the width of the 'comet' tail is smaller 
than dn- For a circular orbit, the large scale appearance is then 
an Archimedes spiral with a step size 27:d„ = crc P orb- 

More generally, the flow speed amplitude and direction 
change along the orbit and the exact shape can differ from a 
perfect spiral. Furthermore, emission along the spiral is not 
constant but depends upon the distance d to the pulsar or, al- 
ternatively, upon the time elapsed since the particles left the 
pulsar The flux received by the observer may also be Doppler 
by factors of a few (de)boosted initially, when the flow speeds 
are high (neglected here). Finally, the orbital orientation will 
also change the appearance of the nebula. For an edge-on sys- 
tem (/ = 90°), the nebula will mimic a microquasar bipolar jet. 

The trajectory of material leaving the pulsar wind termina- 
tion shock can be found as a function of the elapsed time f. 
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Fig. 4. Left: Evolution of the emission along the pulsar nebula. The pulsar nebula is divided into sections at various distances 
Z with the corresponding synchrotron (black) and inverse Compton (grey) spectra shown. The first section shown is for z - \ 
(dz - 0.01). The section dz is then increased: the fifth spectrum is at z » 1.5, the tenth at z ~ 15, the last spectrum shown is 
at z ~ 1000. The total spectrum of the pulsar nebula is shown as a thick black line. The model parameters are the same as in 
Fig- El The thermal emission at ~1 eV is direct light from the stellar companion. Right: Evolution of the particle distribution as 
the shocked pulsar wind moves out from the system. The sections are the same as those shown in the SED. The initial y - power 
law cools rapidly at high energies due to synchrotron radiation. Adiabatic losses dominate at large distances. 



using polar coordinates (r, 0) in the orbital plane. The distance 
d reached along the nebula is 

d^ f vdt (31) 
Jo 

where v is the flow speed calculated in ^5.21 The radius r with 
the coordinates centered on the massive star is given by 

^d^ + 2dd, cos /3 + dj (32) 

where /3 is the angle at which the shocked material flows and 
ds is the star-pulsar separation at f = 0. If Cy is the unit vector 
along the flow direction and e, is the radial unit vector from the 
massive star to the pulsar location at f = 0, then cos/? = Cr.ey. 
The angle 6 of the material is then 

0^a) + }]-ip (33) 

where w is the periastron argument, 77 is the true anomaly of the 
pulsar at f = (angle measured from the star center, 77 = at 
periastron) and tfr is derived from r sin - d sin /3. 

The full pattern follows from having rj vary with time t ac- 
cording to the usual relationships governing orbital motion (i.e. 
using the eccentric anomaly). Note that v, ds and /3 also vary 
with orbital phase. The result is then projected onto the plane 
of the sky, given the orbital inclination /. 

Fig.|5]shows the spiral pattern computed for LS 5039 with 
E = 10^^ erg s"', 7 = 10^, cr = 0.01 and the pulsar located 
at periastron. The figure represents the emission from particles 
'ejected' over a timescale of about an orbital period and a half. 
The spiral 'turns' around and outwards as the pulsar moves 
along the orbit. The lines illustrate the direction at which ma- 
terial flows out from the binary. Spacing between Unes corre- 
spond to equal intervals of time; in other words, lines are more 
spaced out at periastron due to fast orbital motion (before pro- 
jection on the plane of the sky). For very eccentric orbits (such 



as PSR B 1259-63, see 36. 3> . material will preferentially flow 
out in the direction of the major axis on the apastron side (for 
a radial outflow). There will be a preferred emission axis, even 
when the system is seen pole-on. Another effect illustrated by 
Fig.|5lis the bunching of lines due to the projection on the sky 
plane. A binary with a circular orbit but seen at some inclina- 
tion will also have a preferred emission axis. 



Besides the orbital eccentricity and inclination, a third ef- 
fect that will shape the appearance of the cometary nebula is 
the variation of emission with orbital phase. The amplitude can 
change with orbital phase. But even if the spectrum varies lit- 
tle along the orbit, the locus of peak emission changes as the 
elapsed time depends directly on the absolute scale Rs (hence 
phase). If the standoff point is twice as far away at apastron 
than it is at periastron, then the spatial scale on which the radio 
peaks also doubles. In Fig.|5] the intensity of the 5 GHz radio 
emission is represented by the grayscale coloring. Note how 
far out material on the upper left hand side (corresponding to 
apastron) has a higher intensity than material on the right hand 
side (periastron) although the latter was emitted a shorter time 
ago. 



All these effects combine to create radio maps that are 
far from trivial (see §6). Despite the complexity, there is one 
simple prediction: the step of the spiral should be » crcf orb- 
In principle, the radio nebula of y-ray pulsar binaries could 
be used to measure the average magnetization parameter of 
relativistic pulsar winds. The particular cases of LS 5039, 
LSI -H6r303 and PSR B1259-63 are examined in the next sec- 
tion. 
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6. Application to the observed y-ray binaries 

6.1. LS5039 

The wind of the 06.5V companion star in LS 5039 is assumed 
to be radiatively-driven with a velocity and density given by 

i:) ^" ^^^ ^^^^ 

with /3 X I. Spectroscopic observations of LS 5039 constrain 
M„ to * 10"^ Mq yr"' and v^o ~ 2400 km s^' dMcSwain et all 
l2004l) . The wind is assumed isotropic and purely radial. Since 
the orbit of the pulsar takes it to almost a stellar radius of the 
companion at periastron in LS 5039 and LSI -i-61°303, more 
complex wind launching conditions than those assumed above 
could play a role in setting Rs- A significant wind rotation might 
also influence the location of the standoff" point and the flow 
direction. 

The pulsar spindown power was taken to be £ = 
10^^ erg s-\ close to that measured in PSR B 1259-63 ( jO . 
Rs then varies between 2 10" cm (periastron) and 4 10" cm 
(apastron). The two parameters left are the pulsar wind magne- 
tization cr and Lorentz factor jp. Low values of cr > 10""^ yield 
low downstream magnetic fields. Hence, the inverse Compton 
emission dominates with far too little X-rays emitted to match 
the SED. At the other end, values of cr « 0.1 - 1 lead to overall 
emission inefficiency because electrons move quickly out into 
regions of low radiation fields (the downstream flow speed is 
high). Values of cr in the 0.001-0.01 range yield the best results. 
However, with cr x 0.001, the characteristic size of the radio 
emission (Eq. I30> would be too s mall by a factor 10 t o match 
the observed scale of ~ 10 AU jParedes et al. '. 2000. see be- 
low). The best compromise was cr - 0.01. Coincidentally, the 
downstream flow speed saturates at a value (~ 3000 km s"') 
that is then close to the terminal velocity of the stellar wind 
(~ 2000 km s-i). 

The Lorentz factor yp influences the lower cutoff of the par- 
ticle distribution. Lower values enable the inverse Compton 
emission to extend down to the GeV range, as discussed in 
^4.31 With more low energy electrons, the radio emission is 
also more easily accounted for. A value of jp = 10^ was the 
best compromise: at 10^ the GeV emission is unaccounted for, 
at 10"* the radio emission is overpredicted and the downstream 
particle distribution then extends more than four decades in en- 
ergy. The adopted values are cr = 0.01 and jp = 10^. 

6.1.1. The SED 

The evolution of the SED as the material flows away from 
the pulsar is shown in Fig. |2 The high energy SED is very 
much that described in §4, with the spectral breaks in the syn- 
chrotron and inverse Compton spectra at their expected loca- 
tions. Electrons in the comet tail enable the summed emission 
to extend down to low frequencies and also contribute to emis- 
sion in the GeV domain. The comparison to the observations 
shows reasonable but not perfect agreement. Changes in the 
particle distribution, stellar or pulsar wind could reduce the 
discrepancies but have not been attempted here. The purpose 



at this stage is to validate the model without fine-tuning and to 
compare the y-ray binaries under similar assumptions. 

Fig- El shows the variation of the SED between periastron 
and apastron. Apart from the slight lowering of the synchrotron 
break frequency (see ^4.31 . the spectral energy distributions are 
nearly identical. This agrees with the lack of strong orbital vari- 
ability observed in X-rays or radio. However, an orbital modu- 
lation is unavoidable at TeV energies because of yy absorption 
(see below). 

At y-ray energies, the model has plateau emission across 
the EGRET band and up to the HESS regime, as expected in 
§4. However, (1) the GeV flux is underestimated and (2) the 
high energy break is a factor 2-3 too low than what would 
be required to match the hard HESS power-law (~ 2 TeV). 
LS 5039 is in the Galactic Plane so that the EGRET flux could 
be overestimated (i.e. should be interpreted as an upper limit) if 
the diffuse y-ray background subtracted is underestimated or if 
several unresolved sources contribute to the emission. That the 
measurements are non-contemporaneous should not be a prob- 
lem in the present scenario, where variability is expected to be 
small (except for the possible pair cascade, see below). 

The approximate treatment of the emission geometry may 
be responsible for these discrepancies. Adjusting E, cr or the 
standoff point distance can alleviate both problems (Eq. I22> . 
Tweaking the accelerated distribution of particles might help, 
but lacks justification. There could be a direct contribution from 
acceleration processes in the pulsar magnetosphere. Two other 
avenues can also be explored. 

First, the anisotropy of the stellar radiation field might play 
a role. The outgoing energy of an up-scattered photon depends 
on the interaction angle with the higher energies reached at n. 
Even if the distribution of electrons is isotropic, the mean an- 
gle of interaction of the photons detected by the observer will 
depend on the orbital phase, hence there should be an energy 
eff'ect. This is likely to be small at TeV energies as there is not 
much energy dependence on angle in the Klein-Nishina regime 
(^oderskiet al., 2005). 

A second important eff'e ct is yy ab s orptio n 

jBottcher fc De J Eooi lDuhu4 Eoo3: iBednarekt Eoo3) . 
High energy y-rays are emitted in the immediate vicinity of 
the pulsar where they can then interact with stellar photons to 
produce e^e^ pairs. The threshold y-ray energy is ^ 30 GeV. 
The transmitted flux is shown in Fig.|6las grey lines for both 
periastron and apastron. The calculation assumes that y-ray 
emission is isotropic and occurs at the pulsar location (see 
Ipubus 2006 for details). The opacity depends on orbital phase, 
with a maximum close to periastron passage when most of 
the y-ray emission between 0.1-10 TeV should be absorbed. 
An orbital modulation of the VHE flux measured by HESS 
is therefore expected. Note that this is not a discriminating 
feature with the accretion/jet scenario since any y-ray emission 
occurring within x 1 AU of the companion star in this system 
(for instance, in regions close to the black hole) would be 
significantly modulated. 

The average result of absorption will be a hardening of the 
high energy (HESS range) spectrum. Pairs created by the ab- 
sorption of y-rays can subsequently upscatter stellar photons 
to very high energies, thereby initiating a cascade (not taken 
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into account in the calculation shown in Fig.|5). This can sig- 
nificantly increase the emission below 30 GeV as the absorbed 
energy is redistributed below the interaction threshold if the 
cascade is very strong ( Aharonian et al., 2006 ). This wo_uld 
also introduce or bital variability below 30 GeV dPubusiboOfit 
lBednarel3 However, note that high energy pairs above 

ys will radiate predominantly synchrotron emission instead of 
upscattering photons reducing the extent of the cascade 

jAharonian et all 1200611 

Understanding the exact impact of these two effects will re- 
quire complex modelling, including a better representation of 
the shock geometry. The conclusion is that the present fairly 
rustic model is already able to account for the whole SED ad- 
equately, in the manner described in §4 using straightforward 
arguments. 

6.1 .2. Radio maps 

The total radio intensity does not vary much along the orbit but 
the time/distance at which maximum radio emission is reached 
does depend on the orbital phase. The evolution of the radio 
flux with distance to the pulsar was calculated as above for var- 
ious orbital phases and used together with the geometric model 
for the nebula ( ^5.4> to create emission maps of the system. 
Free-free absorption on the stellar wind becomes optically thin 
at ~ 1 AU (see rough estimate in ^2.4> . close to the distance at 
which radio emission starts. Hence, this effect is neglected. 

Fig.0shows the evolution of the 5 GHz emission as a func- 
tion of orbital phase 0. The first map (periastron = 0) shown 
uses the same parameters as those in Fig. |5] except that the 
emission is smoothed by a Gaussian with FWHM of 0.8 mas to 
approximate the observing setup of Paredes et al. (2000). The 
system is seen at an inclination / - 120°, i.e. 30° below the 
orbital plane (/ = 60° would be 30° above the plane, the de- 
generacy from the radial velocity is lifted by the comparison 
of the curvature in the map to the VLBI observations). The 
contours at 50%, 25% and 10% of the maximum surface in- 
tensity can be compared to the 5 GHz VLBI observation. The 
maps at <^ = - 0.2 5 are strikingly similar to that observed by 
IParedes et alJ ll200 0^. with slightly curved emission extending 
a few mas and the inner emission pointing in a slightly different 
direction. The observed very weak counter-jet would be associ- 
ated with emission emitted at the previous periastron passage. 

The model predicts changes in the radio maps with <p 
(Fig.Q that can be verified by further VLBI observations. An 
animated map gives the distinct impression that blobs are al- 
ternatively ejected to the sides of the system, cunningly mim- 
icking a microquasar. Accurate astrometry should show that the 
position of peak radio emission varies by a few mas. The orbital 
motion produces a 'comma' shape that is more pronounced at 
(f> X Q, when particles emitted at apastron start emitting in the 
radio band. Emission occurs principally along the major axis of 
the sky -projected system orbit. This effect shows up as a bunch- 
ing of the lines representing the trajectories in the mid-plane of 

Fig.m 

On progressively larger scales the emission becomes egg- 
shaped along the same axis, with the tip of the 'egg' on the 



apastron side as projected onto the sky. MERLIN and EVN 
radio observations do res olve emission on a larger scale on 
one side than on the other (fP aredesetani2002h . The observed 
emission is more linear than in the model with i = 60°. A 
greater inclination would help since emission becomes exactly 
line-shaped for an edge-on orbit / = 90°. However, the X-rays 
would then be eclipsed when the pulsar moves behind the com- 
panion; an upper limit using the orbit and assuming the X-ray 
emitting region is a sphere of radius /?, gives a; 75°. The ob- 
served extent of radio emission is also larger than expected. On 
an equivalent map to the EVN observation, the 5 GHz intensity 
decreases to 1% of the peak value at 12 mas on one side and 
6 mas on the other instead of 35 mas and 25 mas. The quanti- 
tative agreement is not as good as on scales ~ crcf orb- 

While the map in Fig. shows basically only emission 
within an orbital timescale, maps on scales » ctcPq^^, sum up 
emission from dozens of orbits, representing only a tiny frac- 
tion of the initial energy budget. The physics on such scales is 
undoubtedly more complex than the simple approach taken in 
^5.21 mixing occurs with the surrounding medium; recollima- 
tion shocks re-energise the electrons; shearing introduces tur- 
bulence; magnetic field changes as it is azimuthally stretched; 
shock conditions at earlier times we re different etc. Note that 
the outflow length of -500 AU in IParedes et all ll2002l) co- 
incides roughly with the distance at which the stellar wind 
merges with the ISM. Incidentally, many of these complica- 
tions are also encountered in the propagation of relativistic jets. 

The conclusion is that the pulsar model provides a very 
good quantitative explanation for the radio emission emitted on 
scales of AU in LS 5039 and provides a qualitative framework 
to understand radio emission on larger scales. Radio observa- 
tions have tremendous potential to validate the model. 

6.2. LSI +61° 303 

The pulsar E and yp were taken similar to those of LS 5039 for 
comparison purposes. A slightly higher value of cr = 0.02 was 
adopted here to better match the observed radio scales (§6.2.4; 
note that Table 2 uses cr - 0.01). 

6.2.1. Stellar Wind 

The stellar companion in LSI -i-61°303 is a BOVe type star 
which significantly complicates modelling over the case of 
LS 5039. Ultraviolet and infrared observations of Be stars show 
that their stellar wind has two components: a fast, polar outflow 
with properties typ ical of radiativ ely-driven winds; and a slow, 
equatorial outflow ("Waters! fl986l) . 

A typical polar wind was modelle d using Eq.l34|with M „ - 
IQ-^ Mq yr-i and v^o = 2000 km s"' JWaters et allll988l) . The 
standoff point then varies between 8 10^' cm and 4 10'^ cm. 
The equatorial wind was modelled as 

V,, = vo(R/R^r~^ p,. = po(R/R*r" (35) 

with « «; 3, vq ~ 10 km s~' and po ^ 10"^' g cm"^ 
from the corn parison to infra red observatio ns of LSI -i-61°303 
dWaters et al.[ll988: Marti & ParedesLll995l) . This corresponds 
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Fig. 7. Orbital evolution of the 5 GHz radio emission from the pulsar nebula in LS 5039. The maps were composed by summing 
up the contributions along the spiral nebula, as calculated for the various orbital phases (p (see Fig.|5]for the (p-Q case). The maps 
were then smoothed by a gaussian kernel with a Full-Width at Half Maximum (FWHM) of 0.8 mas (3 10'^ cm at 2.5 kpc). The 
maps share a common color scale to allow direct comparison, with contours at 50%, 25% and 10% of the maximum orbital radio 
emission. The maps at (p-O-Q.25 are very similar to that observed bv .Paredes et aL(.2000.) . 



to an equatorial mass flux about a 100 times larger than the po- 
lar mass flux. The typical opening angle of the Be star equato- 
rial outflow is 15°. Hot observations show the disc extends to 
5-20 /J* and that the terminal velocity is a few hundred km s ' . 
The standoff" point is then between 5 10'*' cm and 2 10'^ cm, 
smaller than the polar wind values because of the higher den- 
sities. The major diff'erence between the two types of stellar 
winds will be at periastron. 



yield better agreement: the downstream magnetic field is then 
lowered, reducing synchrotron losses and enabling stronger ra- 
dio emission at large z- However, the physical extent of the 
associated radio emission would be smaller. A compromise 
would be to have cr change with orbital phase. The distance of 
the shock to pulsar changes by an order-of-magnitude here 
(compared to a factor 2 in LS 5039) and cr could be envisioned 
to decrease as Rs increases. 



6.2.2. The SED 

The resulting SEDs at periastron and apastron for both stellar 
winds are shown in Fig.|8l They agree well with the expecta- 
tions set out in Table |2 most notably the much lower break 
frequencies for the equatorial outflow compared to the polar 
wind due to the smaller R^. 

Orbital variations in X-rays are expected to be greater than 
in LS 5039. Indeed, observations show variations of a factor 
~10 in soft X-rays (0.1-2 keV) and of smaller amplitude at 
higher energies. T he maximum is inferred to be close to peri- 
astron dHarrison et al.. 2000: Casares et al...2005a) . Consulting 
Fig-El this would seem to favor combining a small Rs at perias- 
tron (equatorial wind) and a large at apastron (polar wind). 
At apastron, the pulsar is 15 7?* away from its stellar compan- 
ion, in a region where the Be equatorial disc could already be 
truncated. 

At y-ray energies, the TeV fluxes are consistent with the 
Whipple upper limits. Small Rs for the equatorial flow result 
in high magnetic fields, conversely lowering the y-ray emis- 
sion. Non-intuitively, TeV emission will be higher at apastron 
than at periastron if the pulsar plows through a dense equatorial 
outflow close to its companion and a more tenuous wind fur- 
ther out. With Rs closer to the pulsar, synchrotron losses largely 
dominate at periastron over IC emission. This can be tested by 
observations with VERITAS or MAGIC. As with LS 5039, the 
GeV fluxes are too low to match the EGRET spectrum. An e^e^ 
cascade is unlikely to play a rol e here as the y y opacity is small 
except very close to periastron ( Dubusl EoOd) . 

At low energies, the radio emission is too steep, primarily 
because of the value of cr. A lower value (cr ^ 0.001) would 



6.2.3. Radio outbursts 

LSI +61°303 has regular radio outbursts during which the flux 
increases by an order of magnitude. The outbursts last sev- 
eral days, peaks at an orbital phase ^ 0.2 - 0.7 and repeat 
on the orbital period. There is also an X-ray orbital modula- 
tion (a factor 10 in the ROSAT hwd, less above) that precedes 
the radio peak by a: 0.1- 0.5 (the exact delay varies be- 
tween outbursts). PSR B 1259-63 has similar outbursts at peri- 
astron, a ssociated with the pu lsar crossing the Be equatorial 
outflow (' Melatos et al.L |T995). Indeed, the Be nature of the 
companion is the main feature distinguishing PSR B 1259-63 
and LSI -H6r303 from LS 5039, which has no radio outbursts. 

Assume the pulsar samples the dense equatorial disc around 
periastron, and a more tenuous, fast wind at apastron. Changes 
in Rs, possibly cr, are expected along the orbit; yet these do not 
suffice to explain the outbursts because of the timescales in- 
volved. If such was the case, the ~ 10-day delay in the phasing 
of the X-ray and radio modulations would presumably be asso- 
ciated to the time it takes for electrons to cool from injection. 
The strongest radio emission would occur at apastron, when 
the shock is far from the pulsar and cr is expected to be lower. 
The elapsed time to radio peak is then greater than the orbital 
period (grey line in Fig.|8} and this would smooth out the flux 
modulation. 

The impact of the Be star equatorial outflow is probably 
somewhere else. With a slow wind, the standoff point Rs ~ 
5 ■ 10'*' cm is very small at periastron. Moreover, the equatorial 
wind is very slow so that the Bondi accretion radius oc 1/v^ 
( ^2.4> is actually greater than For the adopted parameters, 
the accretion radius is » 2 ■ 10" cm at periastron. Matter faUs 
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towards the compact object at a rate « 10'^ g s but is not 6.3. PSR B1259-63 
accreted: the pulsar wind is still able to hold off the inflow of 
material However, the situation changes in a significant 

way from that described until now. 

When the relative motions are fast (polar wind case, or high 
velocity isolated pulsar case), the standoff distance is farther 
out than the Bondi capture radius so that the external ram pres- 
sure is mostly unidirectional. The pulsar creates a long comet 
tail directed away from the maximum ram pressure point. At 
periastron, the standoff" distance is within the Bondi capture 
radius. The external pressure from inflowing material is more 
spherically distributed, as in a SNR plerion, so that the shocked 
pulsar wind cannot break open the surrounding medium on a 
short timescale. 

Therefore, a possible picture would be that the shocked 
pulsar wind is embedded in the equatorial disc at periastron 
(plerion-like) and that the geometry is open (comet-like) at 
apastron, where the Bondi radius is smaller than for both 
an equatorial or a polar wind. High energy radiation is not 
affected much as it is emitted close to the pulsar in both 
cases. The change should be mostly in the low frequency emis- 
sion as particles are constantly injected in a bubble slowly 
expanding in the equatorial outffow. The outburst would be 
associated with emission from this bubble of electrons, fol- 
lowing an earlier suggestion bv .Leahv. t2004) . Detailed mod- 
elling is not attempted here but note that adiabatic expansion 
of bubbles has been shown to reproduce the radio outbursts 
llTavlor & Cregorvl Il984l IParedes et all llQQlt lT,eahvl 120041) 



For PSR B 1259 - 63 th e measured E is 8- 10^^ erg s ' 
(Manchest er et all Il995h . The other pulsar parameters were 
taken to be cr - 0.01 and jp - 10^ to compare with 
the other two systems. Emission by containment of a pul- 



by several authors ('Tavani et al., 
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The purpose here is to put the system in the same context as 
LS 5039 and LSI H-61°303, and explore possible large scale 
radio emission. 

Both an equatorial and a polar stellar wind were consid- 
ered, with the same parameters as those given for LSI -i-61°303 
in ^6.2.11 The periastron distance in PSR B 1259-63 is compa- 
rable to the apastron distance in LSI +61°303. Any equatorial 
outffow is then at or close to its terminal velocity and there is 
little difference in Rs with a polar wind (only a factor 3, see 
Table|2}. The equatorial wind is truncated at ~ 15 - 20 7?* and 
is not considered at apastron. Only a polar wind can contain the 
pulsar wind at apastron. 



and, naturally, has also been proposed for the radio outbursts 
of PSR B1259-63 dConnors etalll2002h . 



6.2.4. Radio map 

The radio emission map expected from the simple case of a po- 
lar wind is shown in Fig.|9lfor (p = 0.5 and / = 120°. (Here, 
(f) -Q corresponds to periastron passa ge, to which ^ .2 should 
be added if using the radio ephemeris. lGregorvl 2002'.) T he ori- 
entation of the orbit is illustrated in Fig. 2 of Casares et all 
ll2005ah . The image does not take into account the radio out- 
bursts but should nevertheless be indicative of the overall mor- 
phology to expect. 

As in LS 5039, the nebula has a 'comma' shape. The physi- 
cal scale is larger than in LS 5039 because of the longer orbital 
period. The major axis is almost in the plane of the sky, defin- 
ing a very clear preferred direction for emission. This is consis- 
tent with observat ions of one-sided radio emission on scales of 
10-100 mas ( Massi et all EoOL .2004.) . The mean position an- 
gle will increasingly vary as one goes to higher spatial resolu- 
tions and the impact of the orbital motion becomes clearer. This 
could explain the almost perpendicular directions seen at 3 mas 
and 10 mas resolution, albeit at different epochs (Massi et al.l 
ll993ll200ll) . lMassi et aD(l2004 also find changes in sti'ucture 
about 50 mas away from the core emission on timescale of days 
during the radio outburst. This could be due to a combination 
of rapidly changing emission angles at periastron passage and 
delayed peak radio emission; a proper model for the radio out- 
bursts would be required to verify this possibility. 



6.3.1. The SED 

The spectral energy distribution at periastron and apastron is 
shown in Fig. ^| for the case of a polar wind. The SED at 
periastron for an equatorial wind (not shown) is very close 
to the one shown here, the difference in Rs being small. 
Unsurprisingly, the SED is essentially the same as that for 
LSI H-61°303 at apastron, the orbital separation being the 
same and other parameters varying little. GeV emission from 
PSR B 1259-63 around periastron should be detected by GLAST 
in the future, if only because of this similitude. The TeV slope 
appears in good agreement with the HESS measurement. The 
variations in integrated ffux seen by HESS around periastron 
are beyond this model's present abilities (which predicts none) 
and will probably require taking into account anisotropic IC 
scattering (see Kirket al., 1999), as well as changes in the cir- 
cumstellar density to properly reproduce. There seems to be an 
inffection in the X-ray spectrum between 1-10 keV that might 
correspond to the break observed with RXTE. Very good agree- 
ment can indeed be obtained when using = 10^. 

The SED at apastron is only slightly different. Despite the 
considerably larger orbital separation, the TeV ffux is only low- 
ered by factors of a few. The photon density is indeed much 
lower, but the scale of the ffow (set by R ,) is much larger than at 
periastron. Particles have a longer time available to interact and 
this compensates for the weaker densities. Overall, the radia- 
tion efficiency is the same as at periastron as the same amount 
of energy is being dissipated regardless of orbital phase. The 
model cannot account in its present simplified form for the ob- 
served order-of-magnitude changes in X-ray luminosity. This 
was noted by Murata et al. (200J5) who proposed that lowering 
cr with pulsar distance might help. 

Emission at apastron occurs on large, resolvable scales. The 
grey dashed lines in Fig.^Jshow the amount of the total ffux 
emitted after an elapsed time At > Porh after injection. This ffux 
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is radiated several hundreds of AUs away and may lead to a 
faint arcsecond scale bowshock structure up to IR frequencies. 

6.3.2. Radio outbursts and map 

In PSR B 1259-63, observations of the radio dispersion mea- 
sure i mply the disc is inclined with respect to the orbital plane 
jMelatos et al.[ll99.4 . The radio outbursts seen before and after 
periastron are presumably associated with the two pulsar disc 
crossings (Ball et al., 1999). Working by analogy, smaller sep- 
arations and a coplanar equatorial outflow would explain the 
single outburst per orbit in LSI -i-61°303. As explained above 
for LSI +61°303, the outbursts are unlikely to be due solely to 
a change in The large increase in radio flux might be re- 
lated to the slow expansion of a bubble of relativistic electrons 
formed in the dense flow. Shock-acceleration of electrons from 
the stellar wind itself may also play a role ( B all et allll999t) . 

It is therefore as difficult as with LSI -i-61°303 to accurately 
predict the shape the radio emission might take. Fig. II llshows 
expectations at = 0.1 assuming PSR B1259-63 only inter- 
acts with a polar wind. The spatial scales should be easy to 
resolve since Rs is large. Unfortunately there are no published 
radio maps to compare with. The rapid revolution of the pul- 
sar around periastron produces an arc with two brighter 'knots' 
where pre- and post-periastron electrons are cooling. These two 
knots gradually drift apart. This takes no account of the radio 
outbursts. If these are localised in the equatorial wind, the im- 
pact might be only an increase in the central (unresolved) com- 
ponent, leaving the much fainter diffuse components described 
here mostly unchanged. The bow shock points away from the 
star in almost the same direction for most of the highly eccen- 
tric orbit. 

7. Conclusion 

The pulsar wind scenario offers a tantalizing way to unite 
LS 5039, LSI H-61°303 and PSR B1259-63 into a new-class 
of young, rotation-powered y-ray binaries. At present, nothing 
rules out that the compact object in LS 5039 and LSI -i-6r303 
is a pulsar. The radial velocities allow it and there are no con- 
clusive signs of accretion. Directly unveiling a radio pulse will 
not be possible due to strong free-free absorption; although a 
difficult endeavour, an X-ray or y-ray pulse would prove the 
scenario. 

A similar relativistic pulsar wind with most of its E x 
10^^ erg s"' energy (cr ~ 0.01) in electrons and positrons of 
Lorentz factor ^ 10^ suffices to explain most properties in 
all three objects. The termination shock occurs at the pressure 
balance point. There, particles in the pulsar wind are scattered 
and accelerated into a non-thermal power-law distribution. The 
maximum energy is set by the size of the shock region and ra- 
diative losses. Downstream values in the shocked pulsar wind 
material (density, magnetic field) are set by the MHD jump 
conditions. Particles radiate synchrotron and inverse Compton 
whilst being advected away in a comet tail nebula. At high 
energies, synchrotron dominates over the Klein-Nishina cross- 
section, giving characteristic break frequencies between a ris- 
ing and flat vFy spectrum in hard X-rays, and between a flat 



and a declining vFy spectrum in VHE y-rays. The VHE emis- 
sion occurs close to the pulsar and an orbital modulation of the 
flux in the TeV range is expected in this scenario, due to ab- 
sorption (pair production) on stellar photons. 

The small-scale radio emission is also reasonably repro- 
duced. Radiation shifts to lower wavelengths as particles lose 
their energy to radiation and conditions in the flow change. 
Ultimately, the particles that gave rise to y-ray radiation near 
the star emit radio far from the binary. The flow changes from a 
cometary shape to a spiral-shaped nebula with a step ~ crcf orb- 
Matching with the observed size of radio nebula in LS 5039 
and LSI -i-61°303 implies cr ^ 0.01. On large scales (several 
turns of the spiral), the radio emission is underestimated by the 
model. Changes in the flow conditions on such scales, notably 
velocity and magnetic field evolution, are undoubtedly more 
complex than the simple assumptions made here. Nevertheless, 
the combination of the peculiar geometry, emission at various 
distances and inclination can conspire to create radio maps with 
seemingly bipolar ejections of material as in X-ray binary jets. 
Predictions are made for LS 5039 that are amenable to obser- 
vation with radio interferometry. 

The present model does not provide an entirely satisfac- 
tory view of the GeV y-rays detected by EGRET in LS 5039 
and LSI +61°303. The fluxes are in both cases underproduced. 
The discrepancy might have an observational side, with fluxes 
overestimated, comprising the combined emission from sev- 
eral sources in the wide instrumental PSF. Much more precise 
measurements will be made by GLAST in the future. On the 
theoretical side, inadequacies in the treatment of the emission 
region might be at work, notably anisotropics. The stellar ra- 
diation is anisotropic, the pulsar wind is anisotropic, the ter- 
mination shock is not spherical etc. If synchrotron losses are 
limited, IC cooling of particles will enhance the GeV yield 
and lower the discrepancy. The mismatch between EGRET and 
HE SS spectra of LS 5039 is also difficult for lepton ic jet mod- 
els jPermer & BottcheilEooilParedes et allEoO^. 

The present model does not either satisfyingly explain the 
nature of the radio outbursts in LSI -h61°303 and PSR B1259- 
63. However, the scenario provides a unifying context. As has 
been argued for both sources, the reason must be found in the 
Be nature of the companion stars. The progression of the pul- 
sar in the slow, dense wind could lead to an enclosed geometry, 
rather than an open comet tail geometry. Radio emission might 
also be enhanced by the acceleration of particles from the stel- 
lar wind itself. There is as yet no coherent picture of the radio 
outbursts of LSI H-61°303 in the context of the accretion-driven 
/ relativistic jet scenario. 

Despite these shortcomings, the explanatory scope of the 
model is large. Once the orbital parameters are set, the only 
variables that come into play are E, cr, jp and R^. The SEDs 
come out naturally for reasonable inputs and the comet-tail 
nebula produces radio emission on scales of crPo^^c and be- 
yond. The pulsar scenario therefore seems like a viable alterna- 
tive to the accretion/jet scenario for LS 5039 and LSI -i-61°303. 
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Fig. 1. Observed spectral energy distributions of LS 5039, 
LSI -h61°303 and PSR B 1259-63. Fluxes have been 
transformed to spectral luminosities assuming distances 



of 2.5 kpc rLS 5039. ICasares et alJ |2005b), 2.3 kpc 



s 



(LS I -H61°303.iGregorv et al.lll979ft and 1.5 kpc (PSR B1259- 
63. IJohnston et al.lll994 ." The optical-UV emission peaking 
at ~ 1 eV in the binaries is the light from the stellar compan- 
ion (uncorrected for absorption by the interstellar medium). 
Radio to X-ray flux values in LS 5039 are ta k en from 
(in order of increasin g freque ncy) Ma rti et all (^19981)'); 
IClark et alJ (EoOl): Martocchia et aL (2005.): in LSI -h 61°3() 3 
from IStrickinan et al.1 (Il998l): iMartf & ParedesI (Il995l): 
iMaraschi et alJ lll98ll) : iLeahv et al.' ('l997'): 'Creiner & Rau^ 
1I2OOII): in "pSR B1259-63 from Johnston et al. (2005); 
ICutri et alJ (l2003l) : iMarraco & Orsattil Jl982l) . The RXTE 
spectrum of PSR B 1259-63 obtained during the 2004 pe- 
riastron passage was kindly provided by B. Giebels. The 
upper envelope of radio measurements in LSI +61°303 and 
PSR B1259-63 corresponds to values during radio outbursts. 
The hard X-ray BATSE points in LS 5039 and LSI -i-61°303 
are from [Harmo n et alJ ( l2004l) : H-the OSSE bowtie is from 
Grove et al. ('1995") and the HE y-ray EGRET boxes are from 
Hartman et al. ( 1999) ; EGRE T upper Umits in PSR B 1259-63 



are from lTavani et all ll 19961). The VHE y-r ay HESS measure- 
ments are from Aharon ian et 

ID l2005a'b') and the Whipple 
upper Umits for LSI -i-61°303 fromiFegan et alj (i2005i). 
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Fig. 2. Timescales at the standoff point as a function of the 
electron Lorentz factor y in LS 5039. The radiative timescale 
^rad ~ ^sync + 'ic (^ohd lines) and the escape timescale fesc 
(dashed lines) are shown at periastron (black) and apastron 
(grey). Synchrotron cooling dominates above ys ~ 2 ■ 10^. 
Inverse Compton cooling dominates below, the transition be- 
tween Thomson and Klein-Nishina regimes occurring at 7kn ~ 
6 ■ 1 0'*. 7?, is 2 ■ 1 " cm at periastron (t/ = 0. 1 AU) and 4- 1 " cm 
at apastron (d - 0.2 AU) (see Tabled). 
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Fig. 5. Large scale appearance of the pulsar nebula around 
LS 5039 at periastron. In practice, limited angular resolu- 
tion smoothes out the emission from what is shown here (see 
Fig. 0. The binary is located at the origin. Its size of a few 
tenths of AU is much smaller than the scale of the map. 
Material originating from the standoff point travels outwards 
following the direction of the comet tail at each orbital phase 
(shown by grey lines; the projected outflow path at apastron 
is indicated by a dark line). Because of the orbital motion, 
the emission winds up into a spiral. To illustrate this, dots are 
plotted at the locations reached by material expelled at regu- 
lar time intervals along the orbit. The total elapsed timescale 
represented corresponds to almost an orbital period and a half. 
The greyscale codes the intensity of the 5 GHz radio emission 
(square root scale). At large distances, the step of the spiral 
tends to crcPorb which is ^ 7 AU for the adopted parameters 
(same as those of Figs. l3Pl . The inclination (/ = 120°, corre- 
sponding to a retrograde pulsar orbit, see §6. 1 .2), and the major 
axis orientation {a> = 226°) used here are consistent with the 
measurements of the companion star's rad ial velocity (see orbit 
schematic in Fig. 6 of lCasares et al.l2005lj) . 
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Fig. 6. Spectral energy distribution of LS 5039 at periastron 
(dark full line) and apastron (dark dashed line). The corre- 
sponding grey lines show the fraction received by the ob- 
server after absorption of the VHE y-rays on stellar photons. 
Parameters are identical to those of Figs. I3I4I (but using Rs - 
4-10" cm at apastron). Changes in the overall SED are limited 
and only minor variability is expected except at VHE energies 
where jj absorption strongly modulates the flux on the orbital 
period. 
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Fig. 9. Left: Map of the 5 GHz radio emission from the pul- 
sar nebula in LSI +61°303 at = 0.5 (/ = 120°). is zer o 
at periastron passage (for which ^radio =0.23. 'Gregorvl l2()o3) . 
Right: Corresponding outflow path (as in Fig.|5]but on a lin- 
ear greyscale) before smoothing (the projected outflow path 
at apastron is indicated by a dark line). The stellar wind is 
assumed to be polar (top panel of Fig. |8|i. The smoothing 
Gaussian kernel has a FWHM of 8 mas (2.8 lO'"* cm at 2.3 kpc). 
Contours are at 50%, 25% and 10% of maximum emission. 
There is a preferred direction to the right, towards apastron. 
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Fig. 10. Spectral energy distribution of PSR B 1259-63 at peri- 
astron (solid line) and apastron (dashed line). The pulsar is as- 
sumed to have E - % ■ 10^^ erg s y = 10^ and cr = 0.01. The 
stellar wind is polar (same parameters as LSI +61°303) with a 
standoff' point at = 3 ■ 10'^ cm (periastron) and 4 ■ 10'"^ cm 
(apastron). With an equatorial, slow wind, the standoff point at 
periastron is at - 10'^ cm, not different enough from the 
polar wind case to significantly change the expected SED from 
that shown above. The grey dashed line show the summed late 
time emission (Af > Porb) at apastron, which occurs on scales 
>0.01pc. 



Fig. 8. Spectral energy distribution of LSI +61°303 at perias- 
tron (full black lines) and apastron (dashed black lines). The 
pulsar is assumed to have E - 10^^ erg s"' and y - 10^ as for 
LS 5039, but a slightly higher cr = 0.02. Top: The steflai- wind 
is polar with a standoff point at 7? ^ - 8-10" cm (periastron) and 
4 • 10'^ cm (apastron). Bottom: The stellar wind is equatorial 
with R^ - 5 ■ 10'° cm (periastron) and 2 ■ 10'^ cm (apastron). 
Grey Unes show the fraction of the emission occurring when 
Af > Porb (the periastron line is visible in the top panel but not 
in the bottom panel). 
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Fig, 11. Left: Map of 5 GHz radio emission from PSR B1259- 
63 close to periastron passage (0 - 0.1, i - 30°), assuming 
only a polar-type stellar wind. Right: Corresponding outflow 
path on the same scale before smoothing (the projected outflow 
path at apastron is indicated by a dark line). The smoothing 
gaussian kernel has a FWHM of 0.05" (1.1 10'^ cm at 1.5 kpc). 
The contours are at 50%, 25% and 10% of the maximum inten- 
sity. The tail extends to a few arcseconds. The spatial scale is 
large given the size of the orbit (hence Rs); with an even larger 
eccentricity than LSI -i-61°303, the preferred direction of the 
nebula towards apastron is very clear. 



